Switching off the magnetic exchange coupling by quantum resonances 
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We clarify the role of quantum-well states in magnetic trilayer systems from majority carrier in 
the ferromagnetic and all carriers in the antiferromagnetic configurations. In addition to numerical 
and analytic calculations, heuristic pictures are provided to explain effects of a capping layer and 
side-layer modulation in recent experiments. This immediately offers answers to two unexplained 
subtle findings in experiments and band-structure calculations, individually. Furthermore, it allows 
a more flexible tuning of or even turning off the interlayer exchange coupling. 
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It has been observed [H, 0] more than twenty years ago 
that the interlayer exchange coupling (IEC) of ferromag- 
netic layers in giant magnetoresistance (GMR) system os- 
cillates between ferromagnetism and anti-ferromagnetism 
as a function of spacer width. This feature is ascribed to 
the Ruderman-Kittcl-Kasuya-Yosida (RKKY) oscillation 
[H, |3] and the formation of quantum well states (QWS) 
at Fermi surface [H-Q- These theories, based on the ideal 
system with infinite system widths, are successful at pre- 
dicting the IEC periods. However, the theoretic values 
for their amplitudes always overestimate. Although the 
inevitable roughness on the interfaces was considered as 
the key factor to this disagreement in most cases Q, the 
treatment of approximating the layer thickness as being 
semi-infinite should also play an important role [{J [Toj ) . 
The ab initio calculation [10j in Co/Cu/Co (100) and 
(010) has revealed that the coupling strength is very sen- 
sitive to the variation of side-layer (SL) thickness, which 
redistributes weightings among different RKKY modes. 
The physical origin in (100) case was ascribed to the dis- 
crepancies between the Co density of state in the bulk 
and that in the thin layer, but in (010) the density of 
states is roughly the same and so it was put forward as 
a subtle challenge. 

In addition, the experimentalists 11-13J found that the 
coupling strength displayed an oscillatory feature as the 
SL thickness was varied. It was discovered recently that, 
with an insulating capping layer and by varying SL thick- 
ness, the scattering properties in the SL edge became 
tunable [111]. In contrast to the usage of a metallic cap, 
the oscillation amplitude got stronger and a novel RKKY 
period was deduced[ll|. The quantum resonance, which 
is altered by the scattering paths reflected between the 
bottom of capping layer and the top of spacer, would be 
crucial to answer (1) the origin of disagreements for the 
noble-metal spacerjj] and (2) the creation of new IEC 
period observed in recent experiments 11]. Another mo- 
tivation for studying the finite-size effect for IEC is that 
it may enable us to either enhance or shut off entirely 
the IEC by proper nano-configurations. A known ex- 
ample of such applications is the diluted magnetic semi- 
conductors (DMS), where a room-temperature ferromag- 



netic phase has been obtained by replacing the magnetic 
impurities with the self-organized nanocolumns 14( . The 
quantum resonance in these nanocolumns, which is real- 
ized at properly chosen sizes, enhances the RKKY cou- 
pling between clusters and achieves a higher transition 
temperature for the DMS[l5l|. 

To realize how the boundaries with the capping layer 
and substrate for finite ferromagnetic layers affect the 
IEC, we start with the GMR consisting of a noble- 
metal spacer, in which the minority carrier are con- 
fined in the spacer and their quantum resonance near 
the Fermi surface will determine the position of peaks 
in the IEC oscillation @, [lj|. The Fe/Ag/Fe (001) tri- 
layer is calculated within the single-electron picture with 
different energy barriers in respective layers (l7j. Be- 
sides, we assume that both the substrate and the cap- 
ping layer play the role of infinite potential barriers to 
simulate the recent experimental setup The ex- 



change coupling is defined as the relative energy differ- 
ence J = (AEafm — AEpm) / A with A being the in- 
terface area. The AEfm / AEafm denotes the system 
energy of normal order for the FM/AFM configuration. 
For the case of FM, this energy is calculated by 
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where E(o~, n) is the energy of n-th bound state for car- 
riers with spin orientation a in the confined trilayer, k\\ 
is the transverse wave vector, and T is the temperature. 
Note that Ep/kp is assigned to denote the Fermi en- 
ergy/momentum for the spacer and qp is the Fermi mo- 
mentum for the majority carrier in SL from now on. 

To describe Fe/Ag/Fe (001) system in interface zone 
center at fcgT = 10~ 4 eu, Eq.([T]) and the coupling 
strength J can be estimated by setting Ep = Aev, 1/kp = 
AA and the potential barrier for majority /minority car- 
rier in SL is 5/3et;[rR [3]. In Fig 111 the spacer width 
D is fixed while J is plotted as a function of SL thick- 
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ness, Dpo- Consistent with the experimental findings, the 
oscillation period is dominated by and equals half of the 
Fermi wavelength of majority carrier in the ferromagnetic 
side layer. Similar to the conventional RKKY oscillation, 
J also displays the power-law decay and approaches the 
result of a semi-infinite system as Dp c — > oo. In Fig 121 we 
study the effect of D on J while the SL width is fixed at 
the values which correspond to one of the peaks and bot- 
toms respectively in FigH] Although the extreme values 
of J occur at roughly the sam e sp acer widths, which is 
consistant with experiments 0, Gil], their amplitudes are 
sensitive to the specific choice of SL thickness, with a 
possible enhancement of doubling the semi-infinite value 
at kpD > 14. Overall, the thicker the D, the larger the 
enhancement. 



M 

a. 
W 




or equivalently QWS could be generated periodically as 
a function of spacer width, while the unbound majority 
state is nearly unaffected. The extreme values of J ap- 
pear as the minority carrier form QWS. However, all real- 
istic systems are of finite size; esp. when in nanoscale the 
confinement and formation of QWS in FM/AFM config- 
urations come from not just the minority but all carriers. 
One interesting consequence is that the amplitude of J 
becomes tunable by adjusting the SL thickness, as shown 
in Figl2l Although being unable to derive this, Ref.j!, 0| 
correctly predicted the special spacer widths that could 
render the extreme values. The schematic diagrams in 
FigE] provide heuristic pictures for understanding how 
SL thickness alters the relative positions of quantum res- 
onances for different carriers and makes such a manip- 
ulation possible. The energy barrier height in the side 
layer is set to be zero/infinity for the majority /minority 
carrier for convenience without loss of generality. 
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FIG. 1: (color online) Coupling energies for Fe/Ag/Fe (001) 
trilayer with a fixing D — 5/&f are estimated for different 
SL widths. Both SL are of equal finite width for the circle 
points. In contrast, the square points represents some experi- 
ment conditions where one of the SL is semi- infinite, which we 
mimic by setting D — 200/&F. The solid/dashed vertical lines 
are the theoretic prediction where the peaks/bottom values 
will occur, on which more discussions will follow in Fig(4] 

In the semi-infinite case which contains noble-metal 
spacer, Oretega et al 0, Q noticed that the minority state 
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FIG. 2: (color online) Coupling energies plotted as a function 
of D for Fe/Ag/Fe (001). The solid line denotes the case 
where both SL are set at Af-Dfc = 200 to mimic the semi- 
infinite system. The dashed/dotted lines are for finite systems 
where IcfDf c — 11.5/14 corresponds to the first peak/bottom 
of the solid line in FigfJ] 
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FIG. 3: (color online) Schematic diagrams that exemplify (a) 
the coherent case, which characterizes the simultaneous ap- 
pearance of quantum resonances for carriers in FM and AFM 
configurations and (b) the incoherent case by adding a quarter 
Fermi wavelength to both SL thickness. All carriers in AFM 
no longer sustain QWS due to the violation of boundary con- 
ditions. The big/small arrow denotes the moment /spin ori- 
entation of SL/carrier. 

The SL thickness in FigJ3ta) allows constructive in- 
terferences for the majority carrier in FM (see middle 
plot) and for both carriers in AFM (bottom plot) config- 
urations. Since these resonances are absent in the semi- 
infinite case, they enhance the system energy of both 
configurations by an amount that is roughly double in 
AFM than in FM. When taking their difference, J turns 
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out to be weaker. In Figj3jb), we add a quarter Fermi 
wavelength to both SL thickness. The resonance is ru- 
ined in AFM configuration, but is kept alive in FM. This 
results in an enhanced J as compared to the semi-infinite 
case. Irrespective to whether the capping layer and sub- 
strate are metallic or insulating, the effective reflection 
coefficient at both edges of the spacer can be calculated 
after the inclusion of additional scattering paths that are 
reflected by them 



R P (r c ,D Fc ) = 



r P + r c e 



2iq F D Fc 



1 + rpr c e 2i ^" D ^ ' 



(2) 



where rp and r c denote the reflection coefficient for the 
majority carrier reflected by SL and the capping layer, re- 
spectively. The coefficient Rap for the minority carrier 
obey a similar equation as Eq.Q except the wavevector 
qF in SL becomes complex. The condition for the forma- 
tion of QWS is R p / Ap e 2lkFD = 1 for majority/minority 

carrier in FM, in contrast to RpRApe 2lkFD = 1 for both 
carriers in AFM configuration. These three relations 
allow us to deduce the necessary condition for the co- 
herent case as Rp = Rap as opposed to Rp = —Rap 
for the incoherent one. The connection between the 
hugely tunable J and our heuristic picture of an incoher- 
ent/coherent appearance of QWS is confirmed by com- 
paring Fig|4]with the SL thicknesses that correspond to 
an enhancement/reduction of J in FigJT] 




FIG. 4: (color online) The argument of Rp/Rap is plotted as 
a function of SL thickness. Since \Rp\ = |-Rap| = 1, the zero 
arguments (in circles) denote the case of Rp — Rap which 
allows coherent resonances. While ±7r arguments (squares) 
are for Rp = —Rap that induce incoherent resonances. 

In order to study the enhancement in more detail, we 
start from the asymptotic form of IEC within the free 
electron model and assume that SL are made by the same 
material 



-Epkp 
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By expanding (rp — tap) , the quantum interferences 
for r\ P , rp, and -2rAprp terms are each characterized 



in top, middle, and bottom plots in Fig|3l Quantitatively, 
\tap\ ~ 1 is much bigger than \rp\ in the semi- infinite 
case, as opposed to being exactly equal when finite and 
sustaining QWS. This implies that the maximum en- 
hancement from the incoherent resonance Rp = —Rap 
could reach as high as four times in the semi-infinite 
case. This prediction is verified in FigJSJa). In Fig|2j 
it was observed that the enhancement ratio got bigger 
as D increases. This can also be explained by examin- 
ing the denominator of Eq.(j3]). For finite SL, the con- 
tribution from the majority carrier in FM needs to be 
modified to Rp/(kpD + 2qpD Fc ) 2 , while that from the 
minority remains the same. In contrast, the contribution 
from both carriers in AFM configuration is changed to 
to -R A pRp/(k F D + q F D Fc ) 2 . 

Let's now concentrate on the incoherent case. When 
D 3> 2(qp/kp)D Fc , the contribution from all carriers 
becomes roughly R 2 AP /(kpD) 2 and equals to that of 
the minority. On the other hand, it is much smaller 
when D <C 2(qp/kp)D Fc . So one should adope a thick 
spacer in order to obtain a large enhancement ratio, 
J(D)/J(D <C 2(q F /k F )D Fc ), which depends on 
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k F D 



kpD + 2q F D Fc 



k F D 



kpD + qpD Fc 



(4) 



In the above case, all carriers are totally reflected at both 
edges and satisfy Rp = —Rap- But when the substrate 
is semi-infinite or the top of GMR is either connected 
to a contact or capped by a metal layer, this side of SL 
should be treated as being semi-infinite. It will adopt the 
coefficients rp and rAp for the semi-infinite case, and 
the path factor for additional resonances all equals to 
kpD + qpD Fc . As a result, Eq.Q is modified to 
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k F D 



kpD + q F D Fc 



(5) 



It is obvious that Eq.© is always less than Eq.Q when 
D is large, which explains why the amplitude of the 
square line in Fig[T]is smaller than that of the circle one. 
Both Eqs.(@| and are consistent quantitatively with 
the numerical results in Fig[5][b). The nearly constant 
shift between Eq.([5j) and our calculation is due to the 
fact that the SL thickness we chose only generates ap- 
proximate incoherent resonance in the system with one 
semi- infinite SL. 

A final implication is that we may have solved the puz- 
zle put forward by Halley et al. in Rcf.[Uj. By cap- 
ping Fc/Cr/Fe, they managed to observe a novel period 
for J(D), long predicted and sought after by the band- 
structure calculations As was argued by the authors 
of Ref. , we provide a concrete physical process to ex- 
plain how this additional confinement can lead to revival 
or suppression of the contribution to IEC from a particu- 
lar channel via modifications of the reflection coefficient 
at Fe/Cr. Equations Q and ([5]) have concentrated on 
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FIG. 5: (color online) Panel (a): The peak values of cou- 
pling strength for Fe/Ag/Fe (001) with different arrangements 
of SL thickness: semi-infinite (circle), of the same thickness 
kpDYe = 5.79/7.85 (diamond/square), and with one obeying 
kpDYc = 5.79/7.85 while the other is semi-infinite (inverse tri- 
angle/triangle). The unit Jo is the coupling strength for the 
semi-infinite case with UfD = 14.5. All five lines are fitted by 
a power-law decay with exponent —1.18, —1.38, —1.56, —2.06, 
and —2.41 from up to bottom. Panel (b): Same plot for the 
enhancement ratio with the same thickness /cf-Dfc = 5.79 for 
both SL (diamond) and with one SL set to be semi-infinite 
(inverse triangle) . The solid and dashed lines are plotted from 
Eqs.(gJ and ©. 



the enhancement effect. Let us now derive in more de- 
tail how the suppression comes about. Mainly, one needs 
to create the coherent conditions mentioned in Figj3] (a) . 
Coupling strength of the square line, which represents 
the coherent case, has the highest decaying rate than the 
other four lines in FigJSJa). Its magnitude can be about 
one and two orders smaller than the circle (for the semi- 
infinite case) and the diamond lines (incoherent case) at 
the largest D value in the figure, respectively. This weak 
coupling that results from this nonconventional RKKY 
oscillation, which is generated by Rp = Rap, can be es- 
timated by converting the plus sign for the third term in 
Eq.Q to minus: 



J(D > 2^D Fc ) 



E F D Fe sin(2fc F D 



4tt 2 



D 4 



(6) 



The reason why this decaying power is higher than those 
obtained numerically in Fig|5ta) is that the asymptotic 
form in Eq.Q, on which Eq. (|6j) is based, is only an ap- 
proximation and the higher-order terms are not rigor- 
ously negligible for trilayers with noble-metal spacer |2Cf. 



Although there are many factors that may affect the 
RKKY oscillation in the real Fe/Cr/Fe sample, our calcu- 
lation confirms that the resonances in both SL have the 
capability to drastically change the IEC strength. An- 
other noticeable property is that GMR has multi-period 
RKKY oscillations, which is common in most experi- 
ments, and these modes can be separately tuned by vary- 
ing the SL thickness; i.e., changing their resonant condi- 
tions. Their relative weightings are also affected by the 
SL and capping layer thickness. We are confident that 
the predictions and explanations above for Fe/Ag/Fe can 
be equally applied to Co/Cu/Co. Nevertheless, a quan- 
titatively match will require a more detailed ab initio 
calculation than the pioneering work by Nordstrom et 
aZ.[I3|. For instance, they calculated the density of states 
for a single Co SL and showed that different thickness 
could lead to redistribution of RKKY modes. However, 
they could not locate the origin of this dependence in the 
(010) crystal direction since the density of states they ob- 
tained was not sensitive to varying SL thickness. Based 
on our predictions, we propose to redo the calculations 
for the whole trilayer; namely, include the spacer and 
both SL. 

In conclusions, we find that the strength of IEC for 
GMR with a capping layer can be modulated by as 
much as two orders of magnitude through careful ar- 
rangement of SL thickness. This is made possible by 
creating coherent (to suppress IEC) or incoherent (to en- 
hance) QWS in FM and AFM configurations. This im- 
mediately offers answers to two outstanding puzzles in 
experiments [Tl| and band-structure calculations flfjl] in- 
dividually. We provided heuristic physical pictures and 
both numerical and analytic calculations to support our 
conclusions. Moreover, the discrepancy between theo- 
retic predictions for IEC strength of trilayers and ex- 
perimental findings, which was ascribed to the inter- 
face roughness 0] , also receives an alternative explana- 
tion. Our mechanism carries another potential merit at 
enabling experimentalists to eliminate the intergranular 
exchange coupling in the hard disc drives pi} and ot her 
future spintronic devices. 
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